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ABSTRACT 

A comprehensive, calculated line list of frequencies and transition probabilities for the 
singly deuterated isotopologue of H^", H^D" 1 ", is presented. The line list, called ST1, 
contains over 22 million rotational- vibrational transitions occurring between more than 
33 thousand energy levels; it covers frequencies up to 18500 cm -1 . All energy levels 
with rotational quantum number, J, up to 20 are considered, making the line list use- 
ful for temperatures up to at least 3000 K. About 15% of these levels are fully assigned 
with approximate rotational and vibrational quantum numbers. The list is calculated 
using a previously proposed, high accuracy, ab initio model and consistency checks 
are carried out to test and validate the results. These checks confirm the accuracy of 
the list. A temperature-dependent partition function, valid over a more extended tem- 
perature range than those previously published, and cooling function are presented. 
Temperature-dependent synthetic spectra in the frequency range — 10000 cm are 
also given. 

Key words: molecular spectra, discrete variable representation, H 2 D + , infrared, 
rotation-vibration transitions, line list, partition function, cooling function. 



1 INTRODUCTION 

H2D + is one of the simplest polyatomic quantum systems. 
It consists of two electrons bound to three nuclei (two pro- 
tons and one deuteron) forming a triangular shape at equi- 
librium. The molecule is an asymmetric prolate top with 
three vibrational modes: breathing (v\), bending (y?) and 
asymmetric stretch (^3). All these modes are infrared active. 
H2D + also possesses a permanent electric dipole moment of 
about 0.6 D due to the displacement of the center of charge 
from the center of mass. Hence pure rotational transitions, 
which occur in the far infrared, can occur. This all contrasts 
strongly with the non-deuterated H^" molecular ion which 
has no allowed rotational spectrum and only one infrared 
active vibrational mode. 

In astrophysical environments, H2D + is formed by sev- 
eral reactions; the main one is 



HT + HD — > H2D 



H2 



(1) 



As this and other similar formation reactions are exother- 
mic, the formation of H2D + in the cold interstellar environ- 
ment is favored. Consequently, the abundance of this species 
can be greatly enhanced compared to the underlying D/H 
ratio. As an isotopologue of Hg~, H2D + is a major participant 
in chemical reactions taking place in interstellar medium. In 
particular, it plays a key role in the deuteration processes in 
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these environments |Millar||2003} . Because H3" lacks a per- 
manent dipole moment and hence cannot be detected via ra- 
dio astronomy, its asymmetrical isotopomers, namely H2D + 
and D2H + , are the most promising tracers of the extremely 
cold and highly dense interstellar clouds. These species are 
the last to remain in gaseous state under these extreme con- 
ditions with observable pure rotational spectra ( jDalgarno et\ 
a/|1973 Herbst and Millar 2008 1 . Of these two isotopomers, 
H2D + is the more abundant and easier to observe and hence 
it is the main species to utilize in such astrophysical inves- 
tigations. 



Although the existence of H^D -1 " in the interstellar 
medium (ISM) and astrophysical objects was contemplated 
decades ago ( Dalgarno et alfl973 van Dishoeck et"a/||1992 \ 
with some reported tentative sightin g (jPhillips et qZ||1985 



Pagani et al 1992 Boreiko and Betz 19931, it is only rel- 



atively recently that the molecule was firmly detected in 
the ISM via one of its rotational transition lines ( |Stark ef| 



al 19991. There have been many subsequent astronomical 
studies of H2D + spectra (Caselli et a/||2003| |Vastel et al 



2004) |Ceccarelli et ap004||Stark et qpb04| |Hogerheijde et 



aZ|2006||Harju et aZ|2006| |Cernicharo et afl2007| |Caselli et al 
20081. In particular these spectra have been used to inves- 



tigate the mid-plane of proto-planetary disks ( jRamos et al 
20071, the kinematics of the centers of pre-stellar cores (van 
der Tak et a/|[2005 1 , and the suggestion of possible use as a 
probe for the presence of the hypothetical cloudlets forming 



the baryonic dark matter ( Ceccarelli and Dominik 2006 ) 



2 T. Sochi, J. Tennyson 



Vibrational transitions of H2D 4 " have yet to be observed as- 
tronomically. 

The potential importance of H2D 4 " spectroscopy for cos- 
mology is obvious as it is a primordial molecular species 
which could have a considerable abundance in the early 
universe. In particular it has a potential role in the cool- 
ing of primordial proto-objects (|Dubrovich|1993 Dubrovich 



and Lipovka|1995"||Dubrovich and"F J artridgc 2000; Gal li and 

Z| |2008| ) . Its significance is high- 



Palla||1998] [Schleicher et 



lighted by a number of studies which consider the role of 
H2D + in spectral distortions of Cosmic Microwave Back- 
ground Radiation and whether this can be used to determine 
the H2D + and deuterium abu ndances at different epochs 
( Dubrovich and Lipovka||l995| ). Finally the H 2 D+ 372 GHz 
line has been considered as a probe for the presence of dark 
matter (Ceccarelli and Dominik 2006). 

Due to its fundamental and astrophysical importance, 
H2D + has been the subject of a substantial number of spec- 
troscopic studies. The first successf ul spectroscopic investi- 
gation of H2D 4 " was carried out by Shy et al ( 1981 1 where 



nine rotational- vibrational transitions were measured in the 
infrared region between 1800 and 2000 cm -1 using Doppler- 
tuned fast-ion laser technique, but no specific spectroscopic 
assignments were made. This was followed by other spec- 
troscopic investigations which include the observation and 
identification of the strong and highly-important rotati onal 
lio — In transition line of ortho-H2D + at 372 GHz by Bo 



gey et al (19841 and Warner et al (19841. Amano and Hi- 
rao| ( |2005] ) measured the frequency of the lio — In line of 
H2D + with improved accuracy, as well as the 32i — 3 2 2 line 
at 646.430 GHz. The three fundamental vibrational bands of 



the H2D + ion were observed by Amano and Watson ( 19841, 



Amano (1985) and Foster et al (19861 using laser spec- 



troscopy. Farm'k et al ( 2002 I measured transitions to over- 



tones 2^2 a nd 2^3 and to combinat ion V2 + P3 in jet-cooled 
H2D + ions. Hlavenka et al (20061 measured second over- 
tone transition frequencies of H 2 D + using cavity ringdown 
absorption spectroscopy. Other spectroscopic investigations 



include those of Asvany et al (20071, who detected 20 lines 



using laser induced reaction techniques, and the recent work 
of |Yonezu et al\ ( |2009[ ), who reported precise measurement 
of the transition frequency of the H2D + 2i2 — In line at 
2.363 THz, alongside three more far-infrared lines of H2D + , 
using tunable far-infrared spectrometry. However, none of 
these studies measured absolute line intensities, although 



the work of Farmk et al ( 2002 ) and Asvany et al ( 2007 1 give 



relative intensities. 

Several synthetic, H2D + line lists have been generated 
from ab initio calculations. Prominent examples include the 



line list of Miller et al ( 1989[ ) and another one generated as 
part of the work reported by Asvany et al (2007). Miller et 



aFs list extends to rotational level J = 30 and covers all 
levels up to 5500cm" 1 above the ground state. These lists 
were used in a number of studies for various purposes such 
as spectroscopic assignment of energy levels and transition 
lines from astronomical observations, and for computation 
of a low-temperature partition function ( Sidhu et q/||1992 |. 
The line lists also played an important role in motivating 
and steering the experimental and observational work in 



this field (Miller et al 1989 Asvany et al 20071. Nonethe- 



less, the previous H2D + line lists suffered from limitations 
that include low energy cut-off and the inclusion of a lim- 



ited number of levels especially at high J. These limitations, 
alongside the development of high accuracy ab in itio models 
of H3" and its isotopologues, including H2D + , by Polyansky 
and Tennyson ( |1999| ), provided the motivation to generate a 
more comprehensive and accurate line list. The intention is 
that the new list will both fill the previous gaps and provide 
data of better quality. 

The line list, which we call ST1, consists of 22164810 
transition lines occurring between 33330 rotational- 
vibrational levels. These are all the energy levels with ro- 
tational quantum number J ^ 20 and frequencies below 
18500 cm" 1 . This line list can be s een as a companion to 
the Ht line list of Neale et al ( 1996 1 which has been exten- 



sively used for astrophysical studies; although for reasons 
explained below, the ST1 line list is actually expected to be 
more accurate. 



2 METHOD 

Vibration-rotation calculations were performed using the 



DVR3D code of Tennyson et al (20041. The code calculates, 



among other things, wavefunctions, energy levels, transition 
lines, dipole moments, and transition probabilities. DVR3D 
uses an exact Hamiltonian, within the Born-Oppenheimer 
approximation, and requires potential energy and dipole sur- 
faces to be supplied as input. In general, it is these which 
largely determine the accuracy of the resulting calculations 
( Polyansky et q/|2003 l. 

The vibrational stage of the DVR3D suite requires an 
accurate model for the variation of electronic potential as 
a function of nuclear geometry. Here we use the H^" global 
potential surface of Polyansky et al (2000), which used the 
ultra-high accurate ab initio data of Cencek et al ( 1998 I sup- 



plemented by extra data points; the surface was constrained 
at high energy by the data of |Schinke et al\ ( |1980[). We added 
the H2D + adiabatic correction term fitted by Polyansky and 



Tennyson ( 1999 I to this surface. We also employed Polyan- 



sky & Tennyson's vibrational mass scaling to allow for non- 
adiabatic corrections to the Born-Oppenheimer approxima- 
tion: we used /i H = 1.0075372 u and /irj = 2.0138140 u for 
the vibrational atomic masses, and ^tjj = 1.00727647 u and 

= 2.01355320 u for the rotational atomic masses. The 
accuracy of this model is assessed below. 

Calculations were performed in Jacobi coordinates 
(r\,r2,6), where ri represents the diatom distance (H-H), 
T2 is the separation of the D atom from the center of mass 
of the diatom, and 8 is the angle between r\ and r 2 . Radial 
basis functions of Morse oscillator type were used to model 
r\ (Tennyson and Sutcliffc 19821, while spherical oscillators 
were used to model r 2 (Tennyson and Sutcliffe 1983 1. Follow- 
ing Polyansky and Tennyson (19991, the Morse parameters 
for n were set to r e = 1.71, D e = 0.10 and u e = 0.0108, with 
20 Gauss-Laguerre grid points. Parameters of the spherical 
oscillator functions were set to a = 0.0 and tj e = 0.0075 
with 44 Gauss-Laguerre grid points, following extensive tests 
on convergence of the vibrational band origins. 36 Gauss- 
Legendre grid points were used to represent the angular mo- 
tion. The final vibrational Hamiltonian matrix used was of 
dimension 2000. 

In the rotation stage, the size of the Hamiltonian ma- 
trix, which is a function of J, was set to 1800( J+l) following 
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Table 1. Sample of the ST1 levels file. The columns from left to 
right are for: index of level in file, J, symmetry, index of level in 
block, frequency in cm -1 , vi, i>2, V3, J, K a , K c . We used -2 to 
mark unassigned quantum numbers. 
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18265.61525 


-2 


-2 


-2 


1 


1 
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5 


2490.93374 








1 


2 
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1 


6 
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1 


7 
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2 


2 






Table 2. Sample of the ST1 transitions file. The first two columns 
are for the indices of the two levels in the levels file, while the third 
column is for the A coefficients in s . 



30589 


29553 


7.99E-04 


19648 


18049 


8.37E-03 


8943 


7423 


5.55E-01 


8490 


7981 


2.18E-03 


20620 


22169 


6.91E-04 


17613 


15937 


5.62E-03 


13046 


11400 


1.15E-00 


20639 


20054 


7.26E-03 


14433 


17117 


2.49E-03 


25960 


28074 


1.92E-03 



tests on J = 3 and J = 15. These tests demonstrated that 
choosing sufficiently large values for the rotational Hamil- 
tonian, although computationally expensive, is crucial for 
obtaining reliable results. Our aim was to obtain conver- 
gence to within 0.01 cm -1 for all rotation- vibration levels 
considered. Our tests suggest we achieved this except, pos- 
sibly, for some of the highest lying levels. For these levels 
our basic model, and in particular our corrections to the 
Born-Oppenheimer approximation, are not reliable to this 
accuracy. 

To compute the intensity of the vibration-rotation tran- 
sitions, a dipole moment surface is required. We used the ab 



initio dipole surface of Rohse et al ( 1994 1 to calculate the 
components of the H2D + dipole. In the DIPOLE3 module of 
DVR3D we set the number of Gauss-Legendre quadrature 
points used for evaluating the wavefunctions and dipole sur- 
face to 50; this choice is consistent with the requirement that 
this parameter should be slightly larger than the number of 
DVR points used to calculate the underlying wavefunctions. 

The final ST1 line list consists of two main files: one 
for the levels and the other for the transitions. These two 
files are constructed and formatted according to the method 
and style of the BT2 water line list ( |Barber et aZ|2006[ ). The 
total amount of CPU time spent in producing the ST1 list 
including preparation, convergence tests, and verifying the 
final results is about 8000 hours. We used the serial version of 
the DVR3D suite on PC platforms running Linux (Red Hat) 
operating system. Both 32- and 64-bit machines were used 
in this work although the final data were produced mainly 
on 64-bit machines due to the large memory requirement for 
the high-J calculations. 

This variational nuclear motion procedure used for the 
calculations provides rigorous quantum numbers: J, or- 
tho/para and parity p, but not the standard approximate 
quantum numbers in normal mode, rigid rotor notation. We 
hand labelled those levels for which such quantum numbers 
could be assigned in a fairly straightforward fashion. 5000 of 
these levels are fully designated with 3 rotational (J, K a , K c ) 
and 3 vibrational (vi,V2,V3) quantum numbers, while 341 
levels are identified with rotational quantum numbers only. 
Some of these assignments are made as initial guess and 
hence should be treated with caution. Tables [I] and [2] present 
samples of the ST1 levels and transitions files respectively. 



Table 3. Main experimental data sources used to validate the 
ST1 list. Columns 2 and 3 give the number of data points and 
the frequency range of the experimental data respectively, while 
the last three columns represent the minimum, maximum, and 
standard deviation of discrepancies (i.e. observed minus calcu- 
lated) in cm -1 between the ST1 and the experimental data sets. 




1837 - 1953 
2S39 - 3179 
2839 - 3209 
1837 - 2603 
4271 - 4539 
2946 - 7106 



-0.014 
-0.315 
-0.024 
-0.134 
0.046 
0.008 



0.116 
0.054 
0.054 
0.213 
0.172 
0.242 



0.052 
0.065 
0.019 
0.067 
0.050 
0.088 



3 RESULTS AND VALIDATION 

3.1 Comparison to Experimental and Theoretical 
Data 

We made a number of comparisons between the ST1 line list 
and laboratory data found in the literature. This enabled us 
to validate our results. The main sources of laboratory mea- 
surements that we used in our comparisons are presented in 
Table [3] The table also gives statistical information about 
the discrepancy between our calculated line frequencies and 
their experimental counterparts. Table [4] presents a rather 
detailed account of this comparison for a sample data ex- 
tracted from one of these data sets, specifically that of |As- 



vany et al (20071. This table also contains a comparison of 



relative Einstein B coefficients between theoretical values 
from ST1 and measured values from Asvany et al. The the- 
oretical values in this table are obtained from the calculated 
A coefficients using the relation 



Bi u — 



(2 J' + l)c 3 A ul 
(2 J" + l)i87rhv 3 



(2) 



where A u i and Bi u are the Einstein A and B coefficients 
respectively for transition between upper (u) and lower (I) 
states, J' and J" are the rotational quantum numbers for 
upper and lower states, h is Planck's constant, and v is the 
transition frequency. 

As seen, the ST1 values agree with the measured coef- 
ficients of Asvany et al to within experimental error in all 
cases. Other comparisons to previous theoretical data, such 
as that of [Polyansk y and Tennyson] ( |1999[ ) , were also made 
and the outcome was satisfactory in all cases. 
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Table 4. Comparison between measured ( |Asvany et gf|2007^ and 
calculated (ST1) frequencies and relative Einstein B coefficients 
for a number of transition lines of H2D+ These coefficients are 
normalized to the last line in the table. The absolute B coeffi- 
cients, as obtained from the A coefficients of ST1 list using Equa- 
tion [5] are also shown in column 5 as multiples of 10 14 and in 
units of str.s.kg -1 . 



Transition 




Freq. 


(cm- 1 ) 


B 


Relative 


B 


Vib. 


Rot 




Obs. 


ST1 


ST1 


Obs. 


ST1 


(0,3,0) 


lio «- 


111 


6303.784 


6303.676 


8.05 


0.29 


0.29 


(0,3,0) 


loi <- 


Ooo 


6330.973 


6330.850 


8.59 


0.32±0.02 


0.31 


(0,2,1) 


0oo «- 


111 


6340.688 


6340.456 


7.36 


0.27±0.03 


0.27 


(0,2,1) 


202 <- 


111 


6459.036 


6458.794 


9.17 


0.35±0.04 


0.34 


(0,2,1) 


111 <" 


Ooo 


6466.532 


6466.300 


27.3 


1.00 


1.00 



3.2 Partition Function 

The partition function of a system consisting of an ensemble 
of particles in thermodynamic equilibrium is given by 



i(T) = ^(2Ji + l) gi e 



~ ) 



(3) 



where i is an index running over all energy states of the 
ensemble, Ei is the energy of state i above the ground 
level which has rotational angular momentum Ji , k is Boltz- 
mann's constant, T is temperature, and gi is the nuclear 
spin degeneracy factor which is 1 for para states and 3 for 
ortho ones. 

Using the ST1 energy levels we calculated the partition 
functions of H2D + for a range of temperatures and com- 



pared the results to those obtained by Sidhu et al (19921. 
Table [5] presents these results for a temperature range of 5 
- 4000 K. The results are also graphically presented in Fig- 
ure [l] The table and figure reveal that although our results 
and those of of Sidhu et al agree very well at temperatures 
below 1200 K, they differ significantly at high temperatures 
and the deviation increases as the temperature rises. This 
can be explained by the fact that ST1 contains more energy 
levels which contribute increasingly at high temperatures. 
We therefore expect our partition function to be the more 
reliable one for T > 1200 K. 

Using a Levenberg-Marquardt nonlinear curve-fitting 
routine, we fitted the ST1 partition function curve to a fifth 
order polynomial 



r(T)=^ ai r 



(4) 



in temperature and obtained the following coefficients for T 



in K: 








do 


-0.300315 






'.' = 


+0.094722 






a.2 = 


+0.000571 






a.3 = 


-3.24415 x 


10" 


7 


0,4 = 


+2.01240 x 


10" 


10 


as = 


-1.94176 x 


10" 


14 



(5) 



This polynomial is not shown on Figure [T] because it is 
tually identical to the ST1 curve. 



Table 5. Partition functions of H2D+ for a number of tempera- 
tures as obtained from |Sidhu et q]| | |1992| and ST1 line lists. 



T(K) Partition Function T(K) 
Sidhu ST1 
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2000 
Temperature (K) 



Figure 1. The H2D+ partition functions of ST1 (upper curve) 
and |Sidhu et q/| (jT992) (lower curve). 



3.3 Cooling Function 

ST1 was also used to compute the cooling function of H2D + 
as a function of temperature. The cooling function W which 
quantifies the rate of cooling per molecule is given by 



W(T) = 



' A u i(E u - E l ){2J u + l)g u e 




(6) 



where it and / stand for the upper and lower levels respec- 
tively, and z is the partition function as given by Eq. [3] Fig- 
ure [2] graphically presents our cooling function as a fu nction 
of temperature alongside the H^" cooling function of Neale 
[eFd||T996l ). As seen, the two curves are close for T > 600 K. 
However the H2D + cooling function continues to be signifi- 
cant at lower temperatures, whereas at t hese temperatures 
the cooling curve of was not given by 



Neale et al 



( 1996 1 



because they considered it too small to be of significance. 
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Temperature (K) 



Figure 2. A graph of per molecule emission of H2D+ (upper 
curve) and H^" (lower curve) as a function of temperature on a 
log-log scale. The H2D + curve is obtained from ST1 while the H^ " 
curve is obtained from a digitized image from |Neale et q/| \1996\ . 

This is of course why the cooling properties of HbD" 1 " could 
be important in the early universe. 

Using a Levenberg-Marquardt curve-fitting routine, we 
fitted our HbD" 1 " cooling function curve to a fourth order 
polynomial in temperature over the range 10 - 4000 K and 
obtained 



W(T) = b * Tl 



(7) 



bo = 



-1.34302 x 10" 



fci = +1.56314 x 10" 17 
b 2 = -2.69764 x 10" 19 
6 3 = +7.03602 x 10" 22 

b 4 = -1.10821 x 10" 25 (8) 
The fit was essentially perfect on a linear scale. 

3.4 Synthetic Spectra 

One of the main uses of line lists such as ST1 is to gener- 
ate temperature-dependent synthetic spectra. We generated 
spectra from ST1 using the Spectra-BT2 code, which is de- 
scribed by Barber et q/| ( |2006[> and is a mo dified version of 
the original spectra code of |Tennyson et al\ ( |1993[ ) . Figure [3] 
(a) shows a temperature-dependent H2D + absorbtion spec- 
trum for the region largely associated with pure rotational 
transitions, while Figures [3] (b), (c) and (d) show the cor- 
responding spectrum for the vibrational region. As is usual 
with rotation- vibration spectra, H2D 4 " spectra show a strong 
dependence on temperature. 



4 CONCLUSIONS 

In this study we present a new line list (ST1) for the 
triatomic ion H2D + . The list, which can be obtained from 
the Centre de Donnees astronomiques de Strasbourg (CDS) 
database at [ftp: / / cdsarc.u-stra sbg.fr/pu b/cats /VI / 130[ 



comprises over 33 thousand rotational-vibrational energy 
levels and more than 22 million transition lines archived 
within two files. The tests performed on the line list suggest 
that, although it is based entirely on ab initio quantum 
mechanics, it should be accurate enough for almost all 
astronomical purposes. The one possible exception is for 
predicting the frequency of pure rotational transitions 
which are often needed to high accuracy and which are 
therefore better obtained from measured frequencies using 
combination differences. 
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